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We report a surprising observation that the output directionality from wavelength- scale optical 
microcavities displays extreme sensitivity to deformations of the cavity shape. A variation of the 
cavity boundary on the order of ten thousandth of a wavelength may flip the output directions 
by 180 degrees. Our analysis based on a perturbation theory reveals that a tiny shape variation 
can cause a strong mixing of nearly degenerate cavity resonances with different angular momenta, 
and their interference determines the farfield emission pattern. This work shows the possibility of 
utilizing carefully-designed wavelength- scale microcavities for high-resolution detection and sensing 
applications. 
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Optical microcavities have a wide range of applications 
from lasers, filters, sensors to cavity quantum electrody- 
namics and single-photon emitters pQ. The cavity shape 
has been used as a design parameter to control the out- 
put coupling [2HT4], and a key issue is the sensitivity to 
small variations. In a billiard with closed boundary con- 
dition, the effect of a tiny change of the cavity shape 
accumulates as light undergoes specular reflections in- 
side the cavity. For example, the intracavity ray dynam- 
ics changes from regular to chaotic as the cavity evolves 
from a circle to a stadium shape [15]. However, light in a 
dielectric cavity has a finite lifetime before escaping, and 
the sensitivity to the cavity boundary variation is thus 
reduced. Nevertheless, dramatically different emission 
patterns were observed from similarly deformed cylindri- 
cal polymer lasers [9 j . The observation was explained by 
the distinct geometries of the unstable manifolds, which 
determine the refractive escape routes at large deforma- 
tion. For a smaller deformation from a circle or sphere, 
the output is dominated by tunneling instead of refrac- 
tion [3, 16-20], and the evanescent field becomes highly 
directional even with a weak deformation. Experiments 
using nearly spherical resonators [T6l [21] attributed the 
observed directional tunneling to nonperturbative phase 
space structures in the intracavity ray dynamics. 

All these studies were performed in the semiclassical 
regime, where the cavity size R is much larger than the 
wavelength A. As such, the variation of the boundary, 
though small compared to R, is comparable to or even 
larger than the wavelength. It would be interesting to see 
what happens in the wave regime where R approaches A 
[22U24] . A variation of the boundary with similar size 
relative to R is then on a scale much smaller than the 
wavelength. Can light resolve such a minute structural 
feature, and if so in what way? 

To answer these questions, we preformed numeri- 



cal studies on deformed wavelength-scale dielectric mi- 
crodisks. Our results show that variations of the cav- 
ity boundary on the scale of A/10 4 not only can switch 
the output of a cavity mode from bidirectional to uni- 
directional, but can also flip the output directions by 
180 degrees. To understand such high sensitivity, we ap- 
plied a perturbation theory to the wave problem [25] . 
The shape variation causes a mixing of modes with dif- 
ferent angular momenta, and their interference outside 
the cavity determines the farfield emission pattern. If 
two modes have nearly degenerate frequencies, they can 
be strongly mixed, and the degree of their mixing is ex- 
tremely sensitive to the perturbation of the cavity bound- 
ary. Such extreme sensitivity brings the opportunity of 
utilizing carefully-designed wavelength-scale microcavi- 
ties for high-resolution detection and sensing applica- 
tions. 

In the discussion below we characterize the bound- 
ary of a microdisk cavity by p(0) = R[l + 62 cos(2#) + 
e3cos(30)] in the polar coordinates, where | | , l^l <C 1. 
For a small deformation the dipolar term (ei cos 8) mostly 
leads to a lateral shift of the cavity, and it is eliminated by 
choosing a proper origin of the coordinate system. For a 
wavelength-scale cavity R < A, where A is the wavelength 
in vacuum, or equivalent ly, kR < 2tt where the wave vec- 
tor k = 27r/A. Below we discuss transverse electric (TE) 
modes (electric field parallel to the disk plane) which are 
most common in microdisk lasers, and the same effects 
also exist for transverse magnetic (TM) modes. 

We first consider a series of slightly deformed 
quadrupolar cavities with refractive index n = 3, R = 
1 /im, 62 = —0.01, and varying 63 with | ^3 1 < 10 -3 . Us- 
ing a scattering matrix approach [5j [27] we calculated the 
cavity resonant frequencies and quality (Q) factors. Some 
of the high-Q modes exhibit dramatic output sensitivity 
to the value of 63. Take the resonance at Re[fci2] = 4.387 
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FIG. 1: (Color online) Top row: Internal structure (left), 
farfield pattern (middle), and Husimi Projection on the SOS 
(right) of Mode 1 at kR = 4.387 - il.809 x 10" 5 in a 
quadrupole cavity with R = 1 /mi, €2 = —0.01, £3 = 0, and 
n — 3. The black solid contour represents the farfield ob- 
tained from the 2nd order perturbation theory, which agrees 
almost exactly with the numerical data (red shadow). The 
classical ray dynamics in a billiard of the same shape but 
with closed boundary is also shown on the SOS. The six is- 
lands around sinx = 0.5 correspond to the forward (">") and 
backward ("<") triangular orbits, and the four islands near 
sinx = 0.7 correspond to the diamond orbit ("o"). Bottom 
row: Same as in the top row but the cavity is now sightly 
perturbed with 63 = 10 -4 . The resonance shifts slightly to 
kR = 4.387 - i2.039 x 10" 5 . The classical SOS is mostly 
unchanged, except that the islands of the forward triangular 
orbit disappear. 



most the same. For example, Fig. [T] shows the Husimi 
projection [27j [28] of Mode 1, from which we see that 
it is localized at sinx ~ 0.7 and has little overlap with 
those islands at sinx ~ 0.5. Such high-Q modes can be 
considered as quasi-whispering gallery (WG) modes, and 
their output comes from direct tunneling to the leaky re- 
gion where sinx < 1/n. The tunneling rate is the largest 
where the cavity boundary has the highest curvature, i.e. 
at = 90°, 270° for e 3 = 0. These positions only shift 
by about 1° when e 3 = 10 -4 , which cannot account for 
the strong asymmetry in the farfield pattern. 

To gain insight into the extreme sensitivity of the out- 
put on the boundary deformation, we turn to the angu- 
lar momentum analysis of the modal wavefunction inside 
and outside the cavity, 

' [Z P B p H p (kr) cos(p0), r > p(0). 1 } 

J p (nkr), H p (kr) are the p-th order Bessel function and 
the outgoing Hankel function, respectively. A p (B p ) will 
be referred to below as the Bessel (Hankel) coefficients 
inside (outside) the cavity. Because the cavity has reflec- 
tion symmetry with respect to the horizontal axis, the 
cavity resonances have either even parity or odd parity 
about = 0. Here we consider the even modes, which 
can be decomposed by cos(p0). The analysis of the odd 



(Mode 1) as an example, its output is bidirectional to- 
wards = 0°, 180° at €3 = (see the Supplemental 
Material [25 J; but when e 3 becomes 10 -4 , the emission 
of Mode 1 is greatly suppressed along — 180°, giving 
rise to an enhanced forward (0 = 0°) emitting pattern 
(Fig. [I]). By flipping the sign of e 3 , the output direction 
of Mode 1 is reversed, since now the cavity becomes the 
mirror image of the previous one about the vertical axis, 
i.e. p(tt - 0) = R[l + e 2 cos(2(9) - e 3 cos(30)]. Despite the 
drastic change of the farfield pattern, the spatial pat- 
tern of the resonance inside the cavity barely changes as 
shown in Fig. [I] 

The observed extreme boundary sensitivity cannot be 
explained by semiclassical ray dynamics [3j [4] , which can 
be mapped onto the Poincare Surface of Section (SOS) 
using the positions of rays incident on the boundary 
(represented by the azimuthal angle 0) and the corre- 
sponding angles of incidence x- As shown in Fig. [I] 
the majority of the SOS remains regular in the pres- 
ence of small 62 and e 3 , with unbroken Kolmogorov- 
Arnold-Moser (KAM) curves transversing the entire cav- 
ity boundary (0 G [0, 360°]). There are a few islands cor- 
responding to stable periodic orbits. When e 3 changes 
from to 10 -4 , the islands at sinx °° 0.5 display a no- 
ticeable change. However, the high-Q resonances stay 
at significantly higher sinx, where the SOS remains al- 
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FIG. 2: (Color online) Amplitude of normalized Hankel (a, c) 
and Bessel (b, d) coefficients outside and inside the cavity 
(red crosses connected by solid line) with R = 1 /mi, n = 3, 
62 = —0.01, 63 = 10 -4 (a, b), 10 -3 (c, d). Purple squares and 
black triangles are given by the 1st and 2nd order perturba- 
tion calculation, respectively. Although the dominant angular 
momentum outside has changed from m = 9 to m = 6 as 63 
increases from 10 -4 to 10 -3 , inside the cavity the m — 9 
component is still the dominant one. Inset in (d) shows the 
mode structure inside the cavityat 63 = 10 -3 , which is almost 
identical to that shown in Fig. [ljat £3 = 10 -4 . 
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modes is similar with cos(pO) replaced by sm(p0). Since 
1^2 1 , 1 63 1 <C 1, a quasi- WG mode has a dominant angular 
momentum component m inside the cavity. For Mode 
1, m = 9 [Fig. [2^b, d)]. The quadrupolar deformation 
62 cos(2#) scatters light from m to m±2. Since the m + 2 
component is confined within the cavity more strongly, m 
and m — 2 components are dominant outside the cavity at 
e 3 = 0. They interfere destructively in the = 90°, 270° 
directions, and constructively in the = 0°,180° direc- 
tions, giving rise to the bidirectional output [25]. As 63 
becomes nonzero, the cos (30) deformation introduces ad- 
ditional m±3 components, with m— 3 stronger than ra+3 
outside the cavity. Consequently, the dominant Hankel 
coefficients are m,m — 2,m — 3 as shown in Fig. |5|a) 
at 63 = 10 -4 ; they not only have comparable ampli- 
tudes but also similar phases. Since cos(6#) is symmetric 
about the vertical axis while cos(7#) and cos(9#) are an- 
tisymmetric, it interferes negatively with the other two 
in the = 180° direction while interfering positively in 
the = 0° direction, causing the unidirectional emission 
shown in Fig. [I] 

In the meanwhile, the corresponding Bessel coefficient 
Aq remains much smaller than Ag [Fig. [2^b)] , thus it 
barely alters the mode structure inside the cavity. This 
holds true when 63 increases to 10 -3 , at which Bq domi- 
nates over Bq and Bj outside the cavity. The increasing 
amplitude difference of even and odd Hankel waves re- 
duces the interference effect, bringing down the unidirec- 
tionality, which can be measured by U = dOI(0) cos 
from the normalized farfield intensity 1(0). U is zero for 
isotropic or bi-directional emission, and positive (nega- 
tive) for unidirectional emission in the forward (back- 
ward) direction. Fig. |3ja) shows U as a function of 63. 
As 63 increases, U of Mode 1 rapidly increases to its max- 
imum of 0.39 at 63 ~ 2.7 x 10 ~ 4 before it decreases grad- 
ually. 
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FIG. 3: (a) U versus 63 for Mode 1 (solid line) and 2 (dashed 
line) in Fig. |4ja). 1st order (dotted line) and 2nd order 
(dash-dotted line) perturbation results for Mode 1 are also 
shown, (b) Ratio of Hankel coefficients \Bt/B^\ (solid line) 
and \Bq/Bq\ (dashed line) in Mode 1 as a function of £3- Pa- 
rameters used are the same as the bottom row in Fig. [I] 

The analysis above illustrates that the change of out- 
put directionality is a coherent wave effect. What is 
surprising is the rapid growth of Bq with 63. Fig. [3^b) 
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FIG. 4: (Color online) (a) Complex resonance frequency 
kR (red crosses) in a cavity with R — 1 /xm, 62 — —0.01, 
63 = 10 -4 , and n — 3. The corresponding resonances in a 
circular disk of the same R are marked by black dots. Cyan 
triangles are given by the 2nd order perturbation theory, (b) 
Unidirectionality U of the high-Q resonances (black squares) 
and their distances to the nearest low-Q resonances in the 
complex frequency plane (red triangles) versus Re[kR] of the 
high-Q modes. 

shows the ratio \Bq/Bq\ and \Bj / B$\ as a function of 63. 
B 6 already surpasses B 9 in amplitude at 63 as small as 
2.2 x 10 -4 . How can such a tiny perturbation [es cos(3#)] 
of cavity boundary cause a strong mixing of even and odd 
angular momenta? To identify the cause, we noticed that 
there is a lower-Q resonance whose frequency is close to 
Mode 1. Fig. [4ja) shows a series of high-Q modes with 
nearly constant frequency spacing, and a low-Q series 
with slightly larger spacing. Fig. [4^b) reveals a correla- 
tion between unidirectionality of a high-Q mode and its 
spacing to the nearby low-Q mode. At 63 = 10 -4 , Mode 
1 has the largest U among all the high-Q resonances in 
Fig.|4^a), and its distance to its quasi-degenerate partner 
(Mode 1') is also the shortest. Mode V has a dominant 
angular momentum m' = 6 [25], which appears in Mode 
1 at 63 7^ 0. These results suggests a coupling between 
Mode 1 and V . Recent studies [22] [29] demonstrate that 
a high-Q mode can acquire unidirectional emission from 
a low-Q mode to which it couples. However, this scenario 
does not happen here, because Mode V emits more or less 
symmetrically in the forward and backward directions as 
shown in Ref. [25 . The nature of mode interaction here 
is completely different. 

To gain physical insight about the extreme boundary 
sensitivity, we adopted a perturbation theory [26, 30[l3T] 
to the TE modes in a slightly deformed microdisk cav- 
ity, p(0) = R + e/(0), where the dimensionless e <C 1 
[25] . Since the cavity is nearly circular and the res- 
onances under consideration are quasi- WG modes, we 
chose the resonances fco of a circular cavity of radius 
R as our unperturbed basis and treated the deforma- 
tion 62 cos(2#) + 63 cos(3#) = ef(0) as the perturba- 
tion. fc is determined by the boundary condition for 
TE modes in a circular cavity, which gives T m (koR) = 
(l/n)J^(nk R)/Jm(nk R) - H' m (k Q R) / H m (k^R) = 0. 
In the deformed cavity, the resonance frequency can be 
expanded as fc = fco + fcie + fc2e 2 + 0(e 3 ). For convenience, 
we write A p = a p / J p (nkR), B p = (a p +b p ) / H p (kR) in ([T]) 
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and normalize i/j(r) by scaling the dominant a m to unity. 
In Ref. [25] we show that all a p ^ m and b p are at least of 
order e 1 , thus we define a p ^ m = a p e + /3 p e 2 + 0(e 3 ) and 
6 P = /i p e + 7 P e 2 + 0(e 3 ). By expanding the TE boundary 
conditions [25] to e 2 around r = i?, k = fco, we find the 
corrections to the resonance frequency k as well as the 
coefficients a p and 6 p . 

With the second order corrections f3 p ^ m and j p given 
in Ref. [25], the perturbation theory reproduces the nu- 
merical results nicely (see Figs.[l][4|. In fact, the essence 
of the extreme boundary sensitivity of the output direc- 
tionality is already well captured by the first order cor- 
rections 
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as shown in Fig. [3ja). We have defined 

S p (x) = nJ p (nx)/J p (nx) - H' p (x)/H p (x), F pm = 

c p J Q 27r f(0) cos(p6) cos (mO)dO/2iTR (c p = 2 — 5 P ^) and 
dropped the arguments of the Bessel and Hankel 
functions and their derivatives. We note that the first 
order correction to the resonance, k\ = —ekQF mrni 
vanishes unless f(0) changes the average radius (i.e., 
J f(0)dO 0), thus the second order treatment is needed 
to capture the shift of the resonances (see Fig. [4^a)). 
The main sidebands m ± 2 at €3 = are not affected 
by the variation of 63, since they do not couple to each 
other or m via 63. This holds true both inside and 
outside [Fig. [3^b)] the cavity. 

The presence of another resonance k' R with a domi- 
nant angular momentum m! in close vicinity of koR im- 
plies that T m >(koR) w T m >(k' R) = 0. When this occurs, 
the m! component in ^ m \f;ko) is much enhanced via 
a m f, since T~}(koR) ^> 1. This large prefactor amplifies 
the small boundary variation of cos (m — m')9, especially 
when the m! component is leakier (m! < m) and has a 
stronger influence on the wavefunction outside the cav- 
ity. For example, the unperturbed WG resonance cor- 
responding to Mode 1 is k R = 4.388 - zl.226 x 1(T 5 
with m = 9, and there is a quasi-degenerate reso- 
nance of lower-Q at k' Q R = 4.391 - zl.153 x 1(T 2 with 
m' = 6. The factor \T^} (koR)\ = 7.930 is much larger 
than its typical value in the absence of quasi- degeneracy. 
As a result, a m > increases rapidly with F mm / = 63/2, 
so does B m t with respect to B m . Although a m / ap- 
pears in the Bessel coefficient A m i as well, \B m ' / B m \ 
increases much more rapidly than \A m '/A m \ because 
the scale factor \H m (kR)/H m '(kR)\ is much larger than 
\J m (nkR)/J m >(nkR)\. In Mode 1 the latter is almost 
50 times smaller than the former, which explains the al- 
most identical mode structure inside the cavity while the 
farfield pattern changes dramatically with 63. 

Another important factor for the extreme sensitivity 



is the phase of a m /, which differs from a m (= 1) by ir/2 
as given by With an extra tt/2 relative phase com- 
ing from the the asymptotic form of the Hankel func- 
tion in the farfield, i.e. H p (kr — » 00) ex exp(— ipit/2), 
the m! component interferes constructively with the m 
and m — 2 components in the forward direction and 
destructively in the backward direction. For example, 
at e 3 = 10~ 4 , the farfield intensity 1(0) « [cos(9(9) + 
0.8cos(7(9) +0.5cos(6(9)] 2 . 

We also checked another high Q mode, Mode 2 with 
m — 8 in Fig. [4^ a). Its output directionality U ex- 
hibits a lower sensitivity to the change of 63 compared 
to Mode 1 [Figs. |3ja) andgb)]. It is attributed to the 
larger distance to its nearest low-Q mode with m' = 5, 
which causes less perturbation on Mode 2 (\T^}(koR)\ = 
1.924). However, when U of Mode 1 starts to decrease 
at larger 63, U of Mode 2 keeps increasing [Fig.^a)] and 
reaches a maxima of 0.237 at 63 ~ 1.3 x 10 -3 (not shown). 

In summary, we show that the emission of high-Q 
resonances in wavelength-scale microdisk cavities can 
be changed dramatically when the cavity boundary is 
modified on a scale much smaller than the wavelength 
(~ A/10 4 ). Such extreme sensitivity results from mix- 
ing of quasi-degenerate resonances, and it is expected to 
survive, at least partially, with small boundary rough- 
ness. Although the output directionality is affected by 
the boundary roughness, its sensitivity to the modifica- 
tion of the cavity boundary on top of a given surface 
roughness can still be observed in the numerical simu- 
lation [25], as long as quasi-degenerate modes exist in 
the presence of boundary roughness. Our findings may 
have applications in high-resolution detection and sens- 
ing applications, for example, in detecting and controlling 
thermal vibrations when combined with optomechanics 
techniques [32] . 
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